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ABSTRACT 


The input-output characteristics of an ultrasonic testing system 
used for "stress wave factor" measurements are studied by coupling the 
transmitting and receiving transducers face-to-face, without a specimen 
in between. A summary of some of the fundamentals of digital signal 
processing which are utilized is provided. 

The inputs and outputs are digitized and are processed in a 
microcomputer using digital signal processing techniques. The entire 
ultrasonic test system, including transducers and all electronic 
components, is modeled as a discrete-time linear shift-invariant 
system. A digital bandpass filter is introduced to reduce noise effects 
on the output signal. The output due to a broadband input is deconvolved 
with the input to obtain the unit-sample response and the frequency 
response of the discrete-time system. Then the impulse response and 
frequency response of the continuous-time ultrasonic test system are 
estimated by interpolating the defining points in the unit-sample 
response and frequency response of the discrete-time system. The 
ultrasonic test system is found to behave as a linear-phase bandpass 
filter. 

The unit-sample response and the frequency response of the 
discrete-time model of the test system are used to compute the output 
of the test system to a variety of inputs. Excellent agreement 
between predicted and measured outputs is obtained for rectangular 
pulse inputs of various amplitudes and durations and for tone burst 
inputs whose center frequencies are within the passband of the test 
system. Excellent agreement between predicted and measured outputs 
is also obtained for single cycle inputs of various amplitudes and 
whose center frequencies are within the passband of the test system. 

The input-output limits on the linearity of the system are determined. 


INTRODUCTION 


Conventional ultrasonic testing (UT) is conducted either in the 
through- thickness transmission or pulse-echo modes [1]. In the 
through-thickenss transmission mode, the transmitting and receiving 
transducers are coupled to opposite faces of the structure under 
inspection and the transmitted wave field is analyzed. In the pulse- 
echo mode, the transmitting and receiving transducers, which may 
be combined into a single transducer, are coupled to the same face 
of the structure under inspection and the reflected wave field is 
analyzed . 

Recently Vary et al. [2,3] introduced an ultrasonic nondestructive 
evaluation (NDE) parameter called the "stress wave factor" (SWF) . 
Similar to the pulse-echo test mode, separate transmitting and 
receiving transducers are coupled to the same face of the structure. 
However, unlike conventional pulse-echo testing which is generally 
limited to the analysis of nonoverlapping reflected wave echoes, 
the SWF is also valid for the analysis of overlapping echoes. 
Specifically, an input pulse having a broadband frequency spectrum 
is applied to the transmitting transducer and the number of 
oscillations of the output signal at the receiving transducer exceeding 
a preselected voltage threshold is defined as the SWF. The SWF has 
been correlated with mechanical properties of carbon fiber reinforced 
composites [2-4], Williams et al. [5] theoretically and experimentally 
studied the ultrasonic input-output characteristics of the SWF test 
configuration in a thick isotropic elastic plate. The extension of 
that work to thin plates is currently in progress [6]. 

An important step toward the quantitative analysis of input- 
output relations in any ultrasonic NDE procedure is the quantitative 
characterization of the experimental test system without a test 
specimen. The effects on the output signal due to the ultrasonic 
transducers, the coupling of the transducers to the test specimen, 
and electronic components such as filters, amplifiers, attenuators 
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and cables in the experimental system must be characterized before the 
effects on the output signal due to the test specimen can be isolated. 

This study is part of an overall effort to develop quantitative 
analyses of the SWF and to develop computer-aided nondestructive 
evaluation (CANDE)* capabilities. The experimental UT system is 
characterized by directly coupling the transmitting and receiving 
transducers face to face without a test specimen. Input and output 
signals are digitized using a digital oscilloscope and are processed 
in a microcomputer using digital signal processing techniques. The 
system transfer function of the experimental UT system without any 
test specimen is obtained. 

The results of this study should provide a useful example in the 
characterization of any UT system which utilizes separate trans- 
mitting and receiving transducers. Furthermore, developments in 
CANDE should be facilitated by the digital signal processing procedures 
which are summarized . 


* Pronounced the same as "candy". 
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FUNDAMENTALS OF DIGITAL SIGNAL PROCESSING 


A few results in digital signal processing are summarized in 
this section. Extensive discussions of digital signal processing 
techniques exist elsewhere in the literature [7-10], Primarily, 
the notations in [7] are followed in this outline. 

Digital signal processing is concerned with the representation 
of signals by sequences of numbers and the processing of those 
sequences. Sequences correspond to discrete-time signals derived 
from sampling continuous-time signals. The notation x(n) denotes 
a sequence of numbers whose entries depend on the independent 
parameter n. The sequence x(n) is defined only for integer values 
of n and represents successive samples of a continuous-time 
signal. For example, the sequence x(n) is derived from periodic 
sampling of the continuous -time analog signal x a (t) , where t is time, 
according to 


x(n) = x (nT) (1) 

a 

where T is called the sampling period. The reciprocal of T is called 
the sampling rate or the sampling frequency. The availability of 
high speed digital computers and efficient signal processing algorithms 
have accelerated the implementation of digital signal processing 
techniques compared with continuous-time analog signal processing 
techniques. With the proper analysis, results from the digital 
signal processing of a sequence derived from sampling a continuous- 
time signal can accurately approximate results from continuous- 
time signal analyses. 

A discrete-time system is defined mathematically as a unique 
transformation that maps an input sequence x(n) into an output sequence 
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y(n) . Discrete-time linear shift-invariant systems are discrete-time 
analogs of continuous-time linear time- invariant systems. Table 1 
is a summary of input-output properties of discrete-time linear, 
shift-invariant and linear shift-invariant systems and their 
corresponding analogs in continuous -time systems. Discrete-time 
systems can be further characterized as stable and causal. A stable 
system is one for which every bounded input x(n) produces a bounded 
output y(n). A causal system is one for which the output y(n) for 
n equal to n^ depends only on the input x(n) for n less than or 
equal to n^ . All subsequent discussions here are limited to 
discrete-time linear shift-invariant systems that are stable and 
causal . 

As shown in Table 1, the response of a linear shift-invariant 
system can be characterized by the unit-sample response h(n) via the 
convolution sum given in eqn. (1-1) in Table 1. As a consequence, 
it can also be shown that the steady-state response of a linear 
shift-invariant system to a sinusoidal input is sinusoidal of the 
same frequency as the input but with a magnitude and phase determined 
by the system. It is this property of linear shift -invariant 
systems that makes representations of signals in terms of sinusoids 
or complex exponentials (that is, Fourier representations) so useful 
in linear system theory. Also, it should be mentioned that the 
process of obtaining y(n) from known x(n) and h(n) is called 
convolution, and the process of obtaining h(n) from known x(n) and 
y(n) is called deconvolution. 

Table 2 is a summary of Fourier representations of continuous- 
time signals and of sequences. Fourier representations appear as 
equation pairs whose constituents are often referred to as the 
synthesis and analysis equations as indicated in Table 2. X (jft) 
in eqn. (2-lb) in Table 2 is a continuous function in radian frequency 
ft and is the Fourier transform of the continuous-time signal x (t)» 

cL 

where the symbol j is defined as the square root of minus one. 

If x (t) has units of voltage, X (jft) has units of voltage per 
a a 

radian frequency (or simply volt-sec) . X (jft) is in general complex 
and can be specified by its magnitude and phase. 
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If x (t) is a real-valued function of time, the magnitude and phase 

cL 

of X (ift) will be even and odd functions of frequency, respectively, 
a 

Thus, it is only necessary to present the magnitude and phase of X (jft) 
for positive frequencies. If x (t) Is real and an even function of time, 
the phase of X (j£2) is zero. If x (t) is real and an odd function 

cl cL 

of time, the phase of X (j£2) is ± tt/ 2. If x (t) is shifted (that 

a a 

is, advanced or delayed) in time, the magnitude of X (j^O is not 

cL 

affected but a phase linearly proportional to frequency (that is, 
linear phase) is introduced in the phase of X (jft) . More 

cl 

specifically, the phase (in radians) is proportional to the radian 
frequency with a proportionality constant having units of time. 

The negative of the value of the proportionality constant is the 
time delay [9], 

The remainder of Table 2 deals with the Fourier representations 

of sequences. If properly computed, the Fourier transform X (jQ) 

a 

can be approximated by the Fourier representation of a sequence 
derived from sampling x^(t). Formulae for the Fourier representations 
of sequences are given in the remainder of Table 2. Examples 
of the Fourier representations of continuous-time signals and 
sequences are illustrated in Fig. 1 and will be discussed shortly. 

X(e~ ]a) ) in eqn. (2-2b) in Table 2 is a continuous function of 
radian frequency go and is the (discrete-time) Fourier transform of 
a sequence x(n). The frequency go is in units of radians per 
increment of n, which is simply radians. If x(n) has units of 
voltage, X(e^ W ) has units of voltage per radian. The Fourier 
transform of a sequence is useful for the analysis of general 
sequences. The discrete Fourier series representations of a 
periodic sequence x(n) of period N are also shown in Table 2. 
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A periodic sequence x(n) of period N is one such that x(n) equals 
to x(n + N) . X(k) in eqn. (2-3b) in Table 2 is a periodic 
sequence of period N and the numbers in this sequence are the 
discrete Fourier series coefficients of the periodic sequence 
x(n) having period N. If x(n) has units of voltage, x(k) has 
units of voltage also. The discrete Fourier series of a periodic 
sequence serves as a prelude to the discrete Fourier transform 
(DFT) . X(k) in eqn. (2-4b) in Table 2 is a sequence of length 
N and is called the DFT of a finite-length sequence x(n) of 
length N [7]. If x(n) has units of voltage, X(k) has units of 
voltage also. The DFT is useful for digital signal processing 
because sequences processed by a computer are of finite length and 
also because of the existence of efficient DFT computational 
algorithms . 

By comparing the entries in Table 2 for the discrete Fourier 
series of a periodic sequence of period N and the DFT of a finite- 
length sequence of length N, it is observed that the DFT representation 
is obtained from the discrete Fourier series representation by 
interpreting the finite-length sequence as one period of a periodic 
sequence. The properties of the DFT are similar to those of 
the Fourier transform, except that because of the implied 
periodicity, shifts of x(n) in n by one period N and an integer 
multiple of one period are indistinguishable and a shift in n of 
larger than N is the same as a shorter shift. 

The DFT is computed via a fast Fourier transform (FFT) 
algorithm. To evaluate x(n) from X(k) , the inverse DFT is used 
and is computed via an inverse FFT (IFFT) algorithm. The FFT 
computation is particularly efficient when the length of the 
sequence is an integer power of 2. A sequence of length N is called 
an N-point sequence. The DFT of an N-point sequence is called an 
N-point DFT and is computed via an N-point FFT algorithm. In 
Fig. 1, it will be illustrated that the Fourier transform of a 
continuous— time signal can be reconstructed approximately from 
the DFT of a sequence derived from sampling the continuous-time 
signal. Again, it should be noted that where DFT relations are 
concerned, a finite-length sequence is represented as one period 
of a periodic sequence. 
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Fig. 1 shows Fourier representations of continuous-time signals 
and of sequences derived from sampling the continuous-time signals. 

In part (i) in Fig. la, the continuous-time signal x (t) is assumed 

3 . 

to be bandlimited in frequency with Fourier representation X^Cjfi) 
as shown. The highest frequency of x a (t) is assumed to be Q, q /2. 

In part (ii) in Fig. la, a sequence x(n) results from the periodic 
sampling of x (t) with a sampling period of T. The Fourier transform 

3 

of the sequence x(n) into a continuous function XCe^) is also 

shown. It is observed that X(e Jaj ) can be obtained from X (i£!) by the 

a 

superposition of an infinite number of X^Cjfi) shifted in frequency 
cl) by integer multiples of 2 tt. It is also observed that the magnitude 
and frequency of XCe^) are scaled from those of X (jfi) by 1/T 

3 

and T, respectively. The Fourier transform X (j£!) can be recovered 

3 

exactly, except for the scaling factors, from the Fourier transform 

X(e"^ W ) of the sequence x(n) in the interval of -77 oj <_ tt by 

lowpass filtering XCe-^) using a cutoff frequency of tt, on one 

condition. The condition is the $7 T/2 < ir so that there will be no 

o t 

overlapping of superposed X^(j£2) in forming X(e^ ) . This condition 
is known as the ’’sampling theorem” and assures that if a continuous- 
time signal x (t) is sampled at a frequency greater than twice the 
a 

highest frequency of X (jft) , then X(e JUJ ) is identical to X (jft) in 

3 3 

the interval - tt <_ go <_ tt, except for scaling factors. Because 

X (j£7) is recovered, the continuous— time signal x (t) can be recovered 
a a 

from the sequence x(n) . This mrnimum required sampling frequency 
is called the Nyquist rate or the Nyquist frequency [7, 9, 10]. The 
distortion due to overlapping of superposed X^Cjft) when the sampling 
theorem is violated is called aliasing. In part (ii) in Fig. la, the 
case where there is no overlapping of superposed X (jft) is 

3 

illustrated; thus the signal is not aliased. 

Whereas Fig. la deals with general continuous-time signals. 

Fig. lb deals with finite-duration continuous- time signals. In 
part (i) in Fig. lb, a continuous-time signal x^(t) with nonzero 
values over a finite time duration and its Fourier representation 


7 



X (j£3) are shown. The duration where the signal is nonzero is 
a 

assumed to be t and sampling the continuous-time signal with a 

sampling period T is assumed to result in a sequence of length N. 

It is observed that X a (jf2) for a finite-duration signal 

is not bandlimited in frequency. In part (ii) in Fig. lb, 

a sequence x(n) results from the periodic sampling of x (t) with a 

a 

sampling period of T. All numbers in the sequence x(n) are zero 
for n < 0 and for n > N-l. Thus x(n) is a finite-length sequence of 
length N. The Fourier transform of the sequence x(n) into the 

iti) 

continuous function X(e ) is also shown. Similar to results in Fig. la 

it is observed that X(e~^) can be obtained from X (i^) by the super- 

a 

position of an infinite number of X (j£2) shifted in frequency oo by 

cL 

integer multiples of 2tt. However for this case, because X (jft) is not 

a 

bandlimited in frequency, .there is overlapping of superposed X (jfi) 

• cl 

in forming XCe^) . Thus there is aliasing and X^Cj^) can not be 

recovered exaccly from XCe" 3 ) in the interval of -tt <_ u) _< tt. However, 

if the sampling frequency is greater than twice of most of the 

significant frequency components in X (j£2) , X (jft) can still be 

, a a 

approximated adequately by X(e~* W ) . The DFT X(k) of the finite-length 
sequence x(n) is often computed instead of X(e^ ) . However, in DFT 
evaluations, a finite-length sequence is represented as one period of 
a periodic sequence. 

In part (iii) in Fig. lb, a periodic sequence x(n) of period N 
constructed using the N-point sequence x(n) as a period is shown. 

The discrete Fourier series representation X(k) of the periodic 
sequence x(n) is also shown. X(k) is a periodic sequence of period N. 
The DFT X(k) of the finite-length sequence x(n) can be interpreted as 
a period of X(k) for 0 <_ k <_ N - 1 as shown. It is observed that X(e ) 

be obtained from X(k) by interpolating points in X(k) . Because X (jft) 

i oj . a 

can be approximated adequately by X(e J ) if the sampling frequency 

is greater than twice of most of the significant frequency components 


can 
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in X (ift) , interpolations of X(k) or X(k) can be used to estimate 
a 

X a (jf2) . It is also observed that the numbers in X(k) are scaled 
approximately from the values of X (j£2) by a factor of 1/T and 

cL 

the spacing between successive numbers in the sequence X(k) can 

be interpreted as 2 tt / (NT) in radian frequency which is the 

frequency domain resolution. Again, it is noted that the reconstruction 

of the Fourier transform of a finite-duration continuous -time 

signal from the DFT of a finite-length sequence by interpolation 

is generally an approximation due to aliasing. 

Table 3 is a list of conditions for the special case where 
there is exact correspondence (to within a scaling factor) of the 
Fourier transform of a continuous -time signal and the DFT of a 
finite-length sequence derived from sampling [10] . Also, the N-point 
sequence x(n) can be considered as an M-point sequence if M is 
greater than N with the last M-N numbers in the M-point sequence 
being zero. As discussed here next regarding convolution 
procedures, sometimes it is convenient to extend the length of a 
sequence by appending zeroes. 

Because the response of a linear shift-invariant system can 

be evaluated via the convolution sum in eqn. (1-1) in Table 1, 

convolution representations will be discussed. Table 4 is a summary 

of convolution representations of continuous-time signals and of 

sequences . The Fourier representations of the convolutions are 

also given in Table 4. Eqn. (4-la) in Table 4 gives the linear 

convolution of continuous-time signals x . (t) and x „(t) resulting 

aJ. az 

in x 0 (t). The Fourier representation of the convolution is given 
ai 

in eqn. (4-lb) in Table 4. It is noted that convolution in the 
time domain results in multiplication in the frequency domain. 
Conversely, because of the duality between time and frequency 
domains in Fourier representations, it can be shown that convolution 
in the frequency domain results in multiplication in the time 
domain, such as time windowing [9]. The concepts of periodic and 
circular convolutions will be introduced as preludes to the 
efficient evaluation of linear convolution of sequences using 
DFT procedures. 
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Eqn. (4-2a) in Table 4 gives the periodic convolution of two 
periodic sequences x^(n) and x^(n), each of period N, resulting in 
a periodic sequence x^ (n) of period N. By restricting attention 
to one period of the periodic convolution, circular convolution of 
two N-point sequences x^(n) and x^(n) to obtain an N-point sequence 
x^(n) can be defined from the periodic convolution and is given 
in eqn. (4-3a) in Table 4. The two periodic sequences x^(n) and 
x^Cn), each of period N, are constructed using the N-point sequences 
x-^(n) and x^ (n) as periods, respectively. Eqns . (2-3b) and (2-4b) 
in Table 2 have been used to obtain the simplified eqn. (4-3b) in 
Table 4. This is because x^(n) and x^ 11 ) are equal to x^(n) and 
x^Cn) for n ranging from zero to N-l, thus X^(k) and X^(k) for 
0 < k < N-l become X^(k) and X 2 OO , respectively. The circular 
convolution of N-point sequences is called an N-point circular 
convolution. The designation "circular" is derived from a graphical 
representation of how an N-point sequence can be used to construct 
a periodic sequence of period N [7] . It can be imagined that the 
N points from the N-point sequence are equally spaced in angle around 
a circle with a circumference of exactly N points. The periodic 
sequence of period N is obtained by travelling around the 
circumference of the circle a number at a time recording the N-point 
sequence repeatedly. Also, a rotation of the circle corresponds 
to circular shifting of the N-point sequence [7]. 

The linear convolution of an N^-point sequence x^(n) and an 
N^-point sequence x^ (n) is defined as 

N X"1 

x 3 (n) = l x^(m) x 2 (n-m) (2) 

m=0 


Because the resulting sequence x^Cn) is of length N^+N^-l [7], the 
evaluation of linear convolution via circular convolution requires 
all sequences to be of the same length which has to be N ^+N -1. 
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Eqn. (4-4a) in Table 4 gives the linear convolution of finite-length 
sequences x^(n) and x^Cn) to obtain a finite-length sequence x^(n) 
evaluated via circular convolution. The computation of circular 
convolution using DFT instead of the direct evaluation of linear 
convolution is preferred due to computational efficiencies. The 
(N^+N^-l) -point DFTs of the two convolving sequences are computed 
via a FFT algorithm and then multiplied according to eqn. (4-4b) in 
Table 4. The (N^+N^-l) -point inverse DFT of the product is computed 
via an IFFT algorithm and the result is the desired linear convolution. 
Because the original sequences x^(n) and ^(n) are lengths and 

respectively, N^-l and N^-l zeroes are appended to x^ (n) and X 2 (n), 
respectively, to yield two (N^+N^-l) -point sequences for the (N^+N^-l)- 
point circular convolution. Alternatively, it can be stated that if 
linear convolution of two (N^+N^^l) —point sequences is conducted via an 
(N^+N^-l) -point circular convolution, the last N -1 numbers in one 
sequence to be convolved must be zeroes and the last N.^-1 numbers in 
the other sequence to be convolved must be zeroes. In practice, 
the sequences are extended to a length greater than N^+N^-l to a 
length that Is an integer power of 2 for efficient FFT computation. 

Fig. 2 shows convolution representation of continuous-time signals 

and their sequences derived by sampling. In Fig. 2a, continuous-time 

signals x^(t) of duration and x^Ct) °f duration t^ are convolved 

to form the continuous-time signal x^(t) of duration t^ + t^ . 

In part (i) in Fig. 2b, the continuous-time signal x^(t) is sampled 

with a sampling period T and is assumed to produce an N^-point 

sequence. Then N^-l zeroes are appended to form an (N^+N^-l) -point 

sequence x^(n) . Then a periodic sequence x^(n) of period 

N^+^-l is constructed using x^(n) as a period. Similarly in 

part (ii) in Fig. 2b, the periodic sequence x^Cn) of period 

N,+N 0 -l Is obtained from x 0 (t) which is assumed to produce an N 0 -point 
1 2 a2 2 

sequence when sampled with a sampling period T. Periodic convolution of 
x^(n) and x^Cn) results in the periodic sequence x^(n). Then x^(n) 
can be identified from the circular convolution portion of x^(n) 
by restricting attention to one period of the periodic sequence. It 
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is observed that x^Ct) can be approximated from the interpolation of 
the (N^+N 2 “ 1 )’“P°i nt sequence x^(n). It Is observed that the numbers 
in x^(n) are scaled approximately from the value of x^Ct) by a factor 
of 1/T and the spacing between successive numbers in the sequence 
x^Cn) can be interpreted as T. Because the circular convolution is 
obtained via a DFT, restrictions on the DFT apply to the convolution 
also. Thus the sampling theorem applies to sampling the continuous -time 
signals if x^(t) is to be approximated adequately by the sequence 
x 3 (n) . 

In practice, if a .continuous-time signal is to be sampled, the 
sampling frequency and the sampling length must be specified. Table 
5 is a summary of the significance of the sampling parameters. Eqn. 

C 5—1 ) in Table 5 is written based on the earlier discussion in relation 
to Fig, la. 
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EXPERIMENTAL SYSTEM 


A schematic of the experimental system including the 
specimen is shown in Fig. 3. This is a typical stress wave factor 
test configuration. The goal of this study was to characterize 
the experimental test system without the test specimen. Thus, 
all testing in this study was conducted with the transmitting and 
receiving transducers directly coupled face-to-face without any 
test specimen in between. 

The system consisted of a pulse/function generator (Wavetek 
model 191) ; a 5 MHz (corresponding to -3 dB point) lowpass filter 
(Allen Avionics model F2516) ; broadband (0.1 to 3.0 MHz) transmitting 
and receiving transducers (Panametrics model V105) having an 
approximately flat sensitivity of -85 dB relative to 1 V/yBar; 
an ultrasonic interface couplant (Acoustic Emission Technology 
SC-6); an ultrasonic preamplifier (custom-built by Panametrics) 
having a gain of 40 dB in the frequency bandwidth from 1 kHz 
to 10 MHz; a variable frequency filter (A.P. Circuit Corporation 
model AP 220-5) which could be used as a lowpass, highpass, bandpass, 
or bandstop filter in the frequency range from 10 Hz to 2.5 MHz; a 
digital oscilloscope (Nicolet model 2090 with plug-in model 204-A) 
which could sample and store analog signals at sampling frequencies 
ranging from 0.05 Hz to 20 MHz (corresponding to sampling 
periods ranging from 20 sec to 50 nsec) up to 4096 points; and 
an IBM personal computer (IBM PC) which was interfaced with 
the digital oscilloscope. 

The variable frequency filter was used as a 0.4 MHz to 2.6 
MHz (corresponding to -3 dB points) bandpass filter, selected on 
the basis of the receiving transducer frequency response. 

All signals were sampled at a sampling frequency of 20 MHz 
(that is, one sample every 50 nsec) using 256 points. Thus 
each recorded signal was 12.8 ysec in duration. The 5 MHz 
lowpass filter was selected via antialiasing considerations. 
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Table 6 is a BASIC language program listing for the FFT and IFFT 
computations used in this study. All DFT and inverse DFT evaluations 
in this study were conducted using the FFT and IFFT program in 
Table 6. 
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DIGITAL CHARACTERIZATION OF EXPERIMENTAL SYSTEM 


The objective of this study was to obtain the system transfer 
characteristics between the input to the transmitting transducer 
and the output from the receiving transducer after signal conditioning. 
A schematic of the experimental system is shown in Fig. 3. Once 
again, note that the experiment was conducted with the transmitting 
and receiving transducers coupled face-to-face without a test 
specimen. Thus, the "system" under study consisted of the 
transmitting transducer, the coupling between transmitting and 
receiving transducers, the receiving transducer, the ultrasonic 
preamplifier, the variable frequency filter and all necessary 
cables . 

A continuous-time system can be described by its impulse 
response h (t) and its frequency response H (jft) . In this study, 

H Si 

h (t) and H (jfl!) of the test system were computed from the unit- 
a a 

sample response h(n) and the DFT of h(n) denoted H(k) , respectively, 
of a discrete- time model of the continuous-time system. 

The experimental system was assumed to be a discrete-time linear 
shift-invariant system. As shown in Fig. 3, channel 1 of the 
digital oscilloscope sampled the continuous-time input signal to 
the transmitting transducer. The resulting discrete-time signal 
will be called the input sequence and will be denoted by x(n) . 
Similarly, channel 2 of the digital oscilloscope sampled the 
continuous- time output signal from the receiving transducer after 
preamplification and bandpass filtering. The resulting discrete- 
time signal will be called the output sequence and will be denoted 
by y(n). Thus, the objective of this study was to obtain the 
discrete-time system transfer characteristics between x(n) and 
y(n). The discrete-time system besides being linear shift-invariant 
was observed to be stable and causal. A schematic of the digital 
signal processing procedures applied to obtain the unit-sample 
response and the frequency response of the discrete-time system 
is shown in Fig. 4. 
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The pulse/function generator produced an approximately rectangular 
pulse 20 nsec in duration. The rectangular pulse was passed through 
the 5 MHz lowpass filter before it reached the digital oscilloscope. 
Fig. 5 shows the time histories of the input and output signals by 
linearly interpolating the defining points in the input and output 
sequences x(n) and y(n). As shown in Fig. 5, the durations for 
the input and output signals were approximately 3 ysec and 4 ysec, 
respectively. Thus for the sampling period of 50 nsec, only 
approximately the first 60 and 80 numbers in the sequences x(n) and 
y(n), respectively, were nonzero. Actually, the sequences x(n) 
and y(n) were of length 256 after sampling; that is, the digital 
oscilloscope recorded 256 points each for the input and output signals. 
Then 768 zeroes were appended to the sequences to produce x(n) and 
y(n) of length 1024. The reason for using extended sequences was 
to facilitate evaluation of linear convolution via the DFT as 
discussed later. The DFTs X(k) and Y(k) for x(n) and y(n), 
respectively, were evaluated using a 1024-point FFT algorithm based 
on 1024-point x(n) and y(n) . Because the input and output signals 
were of finite duration, they were not bandlimited in frequency. 
However, because the sampling frequency (20 MHz) was greater than 
most of the significant frequency components in the input signal 
(significant up to 5 MHz) and the output signal (significant up to 
2.6 MHz), the sampling theorem assured an adequate approximation 
of the Fourier transform of the input and output signals by the 
DFT of the input and output sequences. Approximations of the 
magnitude and phase of the Fourier transform of the input and 
output signals are shown in Fig. 5. The magnitude and phase of 
the Fourier transform were obtained by linearly interpolating the 
defining points in the magnitude and phase of the DFT of the 
input and output sequences. The magnitude shown in Fig, 5 is 
presented on a dB scale, normalized with respect to the magnitude 
of the largest component in the Fourier transform. The magnitudes 
of the largest components in the Fourier transforms of the input 
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and output are 2.3 x 10 ^ volt-sec and 9.3 x 10 ^ volt-sec 
corresponding to frequencies of 0.4 MHz and 1.4 MHz, respectively. 
From Fig. 5, it is observed that when the magnitude was small, both 
the magnitude and the phase were erratic because of sensitivity to 
noise contamination. 

Letting the 1024-point DFT of the unit-sample respone h(n) of 
the discrete-time system be defined as H(k) , eqn. (4-4b) in 
Table 4 can then be used to relate the DFT of the input X(k) and 
the DFT of the output Y(k) as 

Y (k) = H(k) X(k) (3) 

where 0 £ k <_ 1023 for this case. Then H(k) for the system can 
be evaluated by rearranging eqn. (3) as 

H(k) = Y(k)/X(k) (4) 

The inverse DFT of H(k) is the desired unit-sample response h(n) . 

The procedure for obtaining h(n) from known x(n) and y(n) is 
called deconvolution. 

However, results from attempts to perform the direct division 
in eqn. (4), using X(k) and Y(k) as shown in Fig. 5, were un- 
satisfactory because of noise. Thus, the output y(n) was filtered 
digitally before division. The digital filter was selected by 
trial and error. Various digital filters were applied and the 
resulting H(k) and the unit-sample response h(n) were used to 
evaluate output properties for various inputs and compared with 
experimental results. The digital filter selected was a bandpass 
filter utilizing time windowing via a 701-point Blackman window [7]. 
Specifically, the digital filter was obtained via the following 

steps: (1) obtain the unit-sample response h , (n) of a desired ideal 

d 

bandpass filter of 0.2 MHz to 4.5 MHz, (2) shift h^(n) by 350 
points to obtain the unit-sample response h^(n), and (3) multiply 
h^(n) by a 701-point Blackman window w^(n) to obtain a 701-point 
unit-sample response h^(n) of the filter. These steps are 
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illustrated in Fig. 4 also. The resulting digital filter was a 

causal linear-phase bandpass filter. Step 1 identified the desired 

bandpass filter. Because an ideal filter is noncausal as indicated 

by its h,(n) which is nonzero for n less than zero, step 2 shifted 
d 

h, (n) by delaying it by 350 points. A causal filter is especially 
important for real-time signal processing. 

Because the desired digital filter was a finite impulse response 

(FIR) filter which by definition should have a unit-sample response 

of finite duration [7], h^(n) was truncated and only numbers between 

n of zero and 701, inclusively, were retained. To minimize the 

effects of truncation, h^(n) was multiplied by a 701-point Blackman 

window w^(n) in step 3 to form h^(n). The length of 701 points 

was found effective in minimizing the effects of truncation and yet 

allowed a 1024-point DFT evaluation. The reason for the length of 

the sequence h^(n) of the filter affecting the DFT evaluation was 

that filtering in the frequency domain (that is, multiplying the 

responses of the filter and the signal in the frequency domain) 

corresponded to convolution in the time domain. Because only 

the first 80 numbers in the sequence y(n) were nonzero, proper evaluation 

of the linear convolution of the sequence y(n) with the sequence 

h, (n) which was nonzero for the first 701 numbers required at lea't 
4 

a 780-point DFT according to discussions associated with eqn. (4-4) 
in Table 4. By increasing the number of points to the nearest integer 
power of 2, which for this case was 1024, 1024-point sequences and 
1024-DFTs resulted. Thus, sufficient zeroes were appended to any 
sequence to extend its length to 1024 points, unless otherwise 
indicated. Also, 1024-point DFTs evaluated via a 1024-point FFT 
were used, unless otherwise indicated. 

(k) denotes the 1024-point DFT of the unit-sample response of 
the digital bandpass filter designed by time windowing the unit- 
sample response of an ideal bandpass filter with a Blackman window. 

Using eqn. (4-4b) in Table 4, the DFT of the output after digital 
filtering was H^(k) Y(k). Then eqn. (4) could be modified as 
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(5) 


H x (k) = H 2 (k) Y(k) / X(k) 

where H^(k) was an approximation to H(k) , and 0 _< k 1023. 

The unit-sample response h^(n) was obtained from H^(k) by 
performing an inverse DFT via an IFFT algorithm. Because h^(n) 
corresponsded to h a (t) which is a real function of time, h^(n) was 
a sequence of real numbers and corresponded to the real part of the 
IFFT. Because a shift of 350 points (corresponding to 17.5 jJsec) 
was introduced in the design of the causal digital bandpass 
filter, the unit-sample response h^(n) had to be shifted back by 
350 points to remove the artifact introduced by the digital filter. 
This was done after the rectangular windowing procedure discussed in 
the next paragraph. 

It was observed that the individual values in the sequence 
h^(n) for n <_ 355 and n 450 were small, less than 1/100-th of the 
maximum value in the sequence. These small values were due to 
noise contamination introduced during signal recording and also by the 
digital filtering. So, the sequence h^(n) was time windowed by a 
rectangular window w 2 (n) defined by 


1 


w 2 (n) 



for 355 < n < 450 
otherwise 


( 6 ) 


Thus, the numbers in h^(n) were not affected for 355 < n < 450, 
but all others had been set to zero. The resulting sequence obtained 
by multiplying h^n) and w 2 (n) is denoted as h 2 (n) . It is observed 
that this time windowing significantly smoothed the phase of H-^k) 
in the high frequency region. 

After the rectangular windowing, the sequence h 2 (n) was shifted 
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back 350 points to remove the shift (that is, time delay) introduced 
by the bandpass digital filter discussed earlier. The 1024-point 
sequence h 2 (n) was shifted to obtain h(n) according to circular 
shifting which in this case becomes 


h(n) = 


1 h 2 (n+350) 

h 2 (n - 1024 + 350) , 
0 


for 0 £ n £ 1023-350 
for 1024-350 £ n £ 1023 
otherwise 


(7) 


The sequence h(n) was the unit-sample response of the test system. 

It is noted that the numbers in h(n) were zero for n > 100 due to the 
windowing described in eqn. (6). Circular shifting instead of 
linear shifting was used because h 2 (n) was the result of a few 
DFT evaluations and when a DFT evaluation was performed, a finite- 
length sequence was implied as a period of a periodic sequence. Thus, 
circular shifting procedures corresponding to shifting a periodic 
sequence were required. 

The DFT of the unit-sample response h(n) was H(k) . Because 
h(n) was nonzero only for < 100, the first 512 points of h(n) were 
used to compute H(k) via a 512-point FFT. Thus, H(k) as obtained 
was a 512-point sequence. H(k) was the frequency response of the 
discrete-time model of the test system. 

Fig. 6a shows the impulse response h (t) of the test system 

3 . 

as obtained by linearly interpolating the defining points in the unit- 
sample response h(n) . An impulse response has units of response 
per unit of excitation multiplied by time [11]. Thus h^(t) in 
Fig. 6a is shown with units of volts per volt-second. It is 
observed that the duration of the unit-sample response was 
approximately 3 Usee. Fig. 6b shows the magnitude and phase of the 
frequency response H^Cjft) of the test system as estimated by linearly 
interpolating the defining points in the frequency response H(k) . 

A frequency response magnitude has units of response per unit of 
excitation [11]. Thus, the magnitude of H^(jfi) has units of volts 
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per volt. The magnitude of H (jfi) is shown in Fig. 6b on a dB scale, 

3 

normalized with respect to the magnitude of the largest component in 

H^(jf2) which is 2.7 x 10 ^ volt/volt corresponding to a frequency 

of 1.7 MHz. It was observed that the normalized magnitude had a 

maximum value at 1.7 MHz and the values of the normalized magnitude 

decreased by 6 dB compared to the maximum value at 0.6 MHz and 

2.3 MHz. It was also observed that the phase of H^Cjfi) was linear 

from 0.3 MHz to 2.7 MHz. The linear phase behavior indicates that 

inputs of frequencies from 0.3 HMz to 2.7 MHz will result in outputs 

delayed in time without distortion. (However, distortion occurs 

because the magnitude of H (jft) is not constant from 0.3 HMz to 

a 

2.7 MHz.) The introduced time delay corresponding to the linear phase 
can be computed to be 0.49 ps. As shown by the magnitude plot, the 
system was insensitive to frequencies below 0.6 MHz and above 2.3 MHz; 
the information in these frequencies was easily contaminated by noise 
and also displayed the effects of digital filtering. Thus, to 
summarize, the system was found to behave as a linear-phase bandpass 
filter in the frequency range of 0.6 MHz to 2.3 MHz. 
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RESULTS AND DISCUSSIONS 


The outputs of the test system due to different inputs can be 
predicted using the unit-sample response h(n) and frequency response 
H(k) of the discrete-time model of the test system. If the input is 
x(n) and the output is y(n) , the input-output relationship is given 
by the convolution in eqn. (1-1) in Table 1. For finite-length 
sequences, eqn. (1-1) in Table 1 can be written in a similar form 
as eqn. (2) . The convolution is conveniently evaluated via the 
DFT. If the DFT of the input is X(k) } then the DFT of the output 
Y(k) can be evaluated using eqn. (3). So, the output y(n) due to the 
input x(n) is the inverse DFT of Y (k) . The frequency response H(k) 
is shown in Fig. 6. Figs. 7-10 show the outputs of the system to various 
inputs. The time histories of the predicted outputs were obtained 
by linearly interpolating the defining points in the sequence y(n) as 
obtained via the DFT procedure using 512-point DFTs . 

Figs. 7 and 8 show the outputs of the system to tone burst 
inputs of approximately 6 psec in duration and 1.5 MHz and 2.5 MHz 
in center frequency, respectively. It is observed that there was 
excellent agreement between the experimentally measured and the pre- 
dicted outputs. The arrival times of peaks of individual cycles 
were predicted to within ± 25 nsec in both cases. The measured peak 
amplitudes of all the individual cycles were correctly predicted 
to within maximum errors of 5 percent and 9 percent for the 1.5 MHz 
and 2.5 MHz inputs, respectively. Similar agreement was also 
obtained for center frequencies ranging from 0.5 MHz to 3 MHz. 

Fig. 9 shows the output of the system to a 7.5 psec duration 
rectangular pulse input. There was excellent agreement between the 
experimentally measured and the predicted outputs. The arrival 
times of peaks of individual cycles were predicted to within ± 25 nsec. 
The measured peak amplitudes of the individual cycles were correctly 
predicted to within a maximum error of 5 percent. From Fig. 9, for 
the approximate time of 5 to 7 psec, both the measured and the 
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predicted output were constant in time. The constant values corresponded 
to 0 and 0.046 volts for the measured and predicted outputs, 
respectively. Similar agreement was also obtained for pulse durations 
ranging down to 0.85 ysec. 

Fig. 10 shows the output of the system to one cycle of 1.75 MHz 
frequency input at 2.3 volts peak-to-peak. The arrival times of peaks 
of individual cycles were predicted to within ± 25 nsec. The shapes 
of the measured and predicted outputs agreed well, but the amplitude 
of the predicted output was larger than the measured output. As 
shown in Fig. 10, the predicted and measured peak-to-peak output 
amplitudes were 8.7 and 7.2 volts, respectively. It was observed that 
the predicted output was larger than the measured output when the 
input and output exceeded 1.4 and 5.1 volts peak-to-peak, respectively. 
This is an effect on nonlinearity and indicates that the linear model 
used in the prediction becomes inadequate for large amplitude inputs 
and outputs. When the input and output were less than 1.4 and 5.1 
volts peak-to-peak, respectively, the measured output peak amplitudes 
were correctly predicted to within a maximum error of 6 percent. 
Restrictions on both input and output amplitudes are imposed because 
the nonlinearity may be due to the input and/or the output. 
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CONCLUSIONS 


The SWF ultrasonic test system input-output characteristics were 
investigated by directly coupling the transmitting and receiving trans- 
ducers face-to-face without a test specimen. Some of the fundamentals 
of digital signal processing were summarized. Input and output signals 
were digitized using a digital oscilloscope, and the digitized data 
were processed in a microcomputer using digital signal processing 
techniques . 

The continuous-time test system was modeled as a discrete-time 
linear shift-invariant system. In estimating the unit-sample response 
and frequency response of the discrete-time system, it was found 
necessary to use digital filtering to remove low amplitude noise which 
interfered with deconvolution calculations. A digital bandpass 
filter constructed with the assistance of a Blackman window and a 
rectangular time window were used. Approximations of the impulse 
response and the frequency response of the continuous-time test system 
were obtained by linearly interpolating the defining points of the 
unit-sample response and the frequency response of the discrete-time 
system. It was found that the test system behaved as a linear-phase 
bandpass filter in the frequency range of 0.6 MHz to 2.3 MHz. 

These frequencies were selected in accordance with the criterion that 
they were 6 dB below the maximum peak of the amplitude of the 
frequency response. 

Furthermore, using the unit-sample response and the frequency 
response of the discrete-time system and the linear shift-invariant 
system assumption, the output of the system to various inputs was 
predicted and the results were compared with the corresponding 
measurements on the system. The predicted output was obtained by 
linearly interpolating the defining points of the discrete-time 
output sequence. Tone bursts of various center frequencies and 
durations, rectangular pulses of various durations, and angle cycle 
inputs at 1.75 MHz with various amplitudes were considered. 
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Excellent agreement between the predicted and measured results was 
obtained for tone bursts with center frequencies from 0.5 MHz to 
3 MHz and rectangular pulses. Excellent agreement between the 
predicted and measured results was also obtained for one cycle inputs 
of various amplitudes , as long as the input and output were less than 
5.1 volts peak-to-peak, respectively, at a 40 dB system gain. 

These results are specific for the particular set of components in the 
test system. 

The results of this study will be useful in the quantitative 
analysis of the SWF when test specimens are inserted between the 
transducers. With known test system characteristics, the effects 
on the output signal due to the test specimen having a variety of 
flaw states can be isolated. Furthermore, by using and discussing 
digital signal processing methods extensively, developments in 
CANDE may be facilitated. 
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TABLE 1 Input-Output Properties of Discrete-Time Linear, 

Shift-Invariant and Linear Shift-Invariant Systems 
and Their Corresponding Analogs in Continuous -Time 
Systems . 


Discrete-Time 

System 

System Input-Output 
Properties 

Analog in 
Continuous-Time 
System 

Linear 

If x^n) -* y^n) 
and x 2 (n) ■+ y 2 (n) , 
then ax^(n) + bx 2 (n) 

■+ ay 1 (n) + by 2 (n) 
(Principle of Superposition) 

Linear 

Shif t-Invariant 

If x(n) ■> y (n) , 

then x(n - k) -*■ y(n - k) 

Time-Invariant 

Linear 

Shift-Invariant 

If 6 (n) -*• h(n) , 
then x(n) -*■ y(n) 

oo 

where y(n) - Z x(k)h(n - k) 

k=J ° * 
(1-1) 

(Convolution Sum) 

Linear 

Time- Invariant 
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TABLE 2 Fourier Representations of Continuous^-Time Signals 
and of Sequences. 


1 

Fourier Representation 

Synthesis Equation 

Analysis Equation 

Fourier transform of 

continuous -time 

signal x (t) 
a 

f 00 

x (t) = X (jJOe^dfi 

a 2 tt J a 

(2-la) 

- CO 

X (jft) = x (t)e ^ fc dt 

a J a 

(2- lb) 

Fourier transform of 
sequence x(n) into 
the continuous 
function XCeJ^) 

x(n) = ■— X(e^ aJ )e^ aJn dw 

-tt (2-2a) 

OO 

X(e^ W ) = J x(n)e 
n=-°° 

(2-2b) 

Discrete Fourier 
series of periodic 
sequence x(n) of 
period N 

x(n) = i Y X(k)e jk(2,r/H)n 
k=0 

(2-3a) 

X(k) = l 5(n)e“ jk(2 ^ /N)n 

n=0 

( 2— 3b ) 

Discrete Fourier 
transform (DFT) 
of finite-length 
sequence x(n) of 
length N into the 
sequence X(k) of 
length N 

x(n) Y X(k)e Jk(2,r/N)n , 

k=0 

for 0 < n < N-l; 
x(n) = 0 , otherwise. 

(2-4a) 

X(k) = Y x(n)e-J k(2ll/N)n 

n= o 

for 0 <_ k <_ N-l; j 

j 

X(k) = 0 , otherwise. ; 

(2-4b) ; 
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TABLE 3 Conditions for the Special Case Where There is 

Exact Correspondence (to within a Scaling Factor) 
of Fourier Transform of a Continuous -Time Signal 
and DFT of a Sequence Derived from Sampling [10] . 


Conditions 

(1) The continuous-time signal x (t) must be periodic. 

a 

(2) The frequency content of x (t) must be bandlimited. 

cL 

(3) The sampling frequency must be greater than twice 

the largest frequency component of x (t) . 

a 

(4) The sampling length must truncate x (t) with a 

a 

truncation interval exactly equal to one period 

(or integer multiple of one period) of x (t) . 

a 


TABLE 4___ Convolution Representations of Continuous-Time Signals 
and of Sequences. 


1 

Convolution Representation 

Definition 

Fourier Representation 
of Convolution 

i 

I 

1 Linear convolution of 

continuous-time signals 

| x al <C> ^ x a2 (t) 

1 

x a3 (t) " | x al CT)x a2 (t - T)dT 

(4-la) 

x a3 (jft) - x al (jfi)x a2 (jfi) 
(4-lb) 

Periodic convolution of 
periodic sequences 
x^(n) and x 2 (n), each 
of period N 

N-l 

x 3 (n) = l x 1 (m)x 2 (n-m) 
m=0 

(4-2a) 

X 3 (k) = X 1 (k)X 2 (k) 

(4-2b) 

Circular convolution of 
finite-length sequences 
x^(n) and X2(n) , each of 
length N. The sequences 
are represented as 
periods of periodic 
sequences x^(n) and 
j^CnJjeach of period N. 

N -i 

x (n) = l x (m)x (n-m) 
m=0 

for 0 _< n <_ N-l; 
x^(n) = 0 , otherwise. 

(4-3a) 


Linear convolution via 
circular convolution of 
finite-length sequences 
x^(n) with length 
and X2(n) with length N2. 
The length of each 
sequence is increased to 
Nj_+N 2-1 by appending 
sufficient number of 
zeroes to each sequence. 
The resulting sequences 
are represented as 
periods of periodic 
sequences x^(n) and 
^(n), each of period 

Ni + N 2 “l 

N x +N 2 -2 

X3 (n) = l x 1 (m)x 2 (n-m) 

m=0 

for 0 < n < N^+N 2 ~2; 
x^Cn) = 0 , otherwise. 

(4-4a) 

X 3 (k) = X x (k)X 2 (k) 

for 0 k ^ N 1 +N 2 -2; 

X 3 (k) = 0 , otherwise. 

(4-4b) 
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TABLE 5 


Significance of Sampling Parameters for Deriving 
a Sequence from Sampling a Continuous -Time Signal. 


Sampling Parameter 

Significance 

Relationship 

Sampling frequency 1/T, 
where T is sampling period 
in sec and 1/T is frequency 
in Hz . 

To avoid distortion 
due to aliasing, the 
sampling frequency 
must be larger than 
twice the highest 
frequency component 
contained in the con- 
tinuous-time signal. 

ft 

— < — (5-1) 

2 T k ' 

where ft 0 /2 is the 
maximum radian 
frequency of the 
c on t inuo us - 1 ime 
signal . 

Sampling length NT, 
where N is length of the 
sequence derived from 
sampling and T is 
sampling period in sec. 

To improve DFT fre- 
quency domain 
resolution, the 
sampling length must 
be increased. 

Frequency 
resolution (Hz) 

= 1/ (NT) 

(5-2) 


To avoid distortion 
due to "leakage" 
[23,25] resulting 
from truncating the 
signal, the sampling 
length must be 
sufficient to totally 
contain the signal?* ** 

NT _> Time duration 
of signal 

(5-3) 


To utilize the best 
efficiency of FFT 
computation, N should 
be an integer power 
of 2. 

N = 2 V (5-4) 

where v is an integer. 


* If the sequence is to encounter a linear convolution evaluation using 
circular convolution via a FFT algorithm, appended zeroes may be 
required to extend the length of the sequence. 


** If the signal is periodic, the sampling length should be equal to 
the period or an integer multiple of the period. If truncation 
can not be avoided, the truncation effects should be minimized 
by smoothing the signal by proper time windowing [7-10]. 
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TABLE 6 


Program Listing for Microcomputer-Based N-Point FFT and 
IFFT Computations in BASIC when N is an Integer Multiple 
of 2. 


30000 
30010 
30020 
30030 
30040 
30050 
30060 
30070 
30080 
30090 
30100 
30110 
30120 
30130 
30140 
30150 
30160 
30170 
30190 
30200 
30210 
30220 
30230 
30240 
30250 
30260 
30270 
30280 
30290 
30300 
30310 
30320 
30330 
30340 
30350 
30360 
30370 
3 0380 
30390 
30400 
30410 
30420 
30430 
30440 
30450 
30460 
30470 
30480 


rem Subroutine to calculate Fast Fourier Transform (FFT) 

rem or Inverse Fast Fourier Transform (IFFT) 

rem 

rem Data are in fr(real) and fi ( imaginary ) arrays 

rem Calculation is in place, output replaces input 

rem The following parameters must be defined : 

rem n££t% : The number of points (must be power of 2) 

rem itr% : 1 for FFT, -1 for IFFT 

rem 

rem Check if nfft# is integer power of 2 
xnf f t=nf f t% 
for i%=1 to 20 

xnf f t= . 5*xnf f t : if xnfft=1 then 30170 
if xnfft>1 then 30160 

print " ERROR in subroutine FFT ( nff t% wrong ! ! ! ) " : print 

stop 

next 

rem Check itr% 

if itr#=1 or itr#=-1 then 30220 

print " ERROR in subroutine FFT ( itr% wrong ! ! ! ) " : print 

stop 

rem Calculate FFT 
mr%=0 : nn%=n££t%-1 
for m#=1 to nn# 
nl#=nf f t% 

nl%=nl%\2 : if mr#+nl#>nn$ then 30260 

mr%= (mr% mod nl%) +nl% : if mr#<=m# then 30300 

tr=fr(mj6+l) : f r ( m#+1 ) =f r ( mr^+1 ) : fr(mr% + l) = tr 

ti=fi(m#+l) : f i ( m%+1 )= f i ( mr^+1 ) : fi(mr$+l)=ti 

next m% 

nl%=1 

if nl%>=n££t% then 30440 
is tep%=nl%+nl% : el=nl# 
for m% = 1 to nl% 

arg=3 . 1 41 593*( itr#*( 1 -m%) ) /el : wr-cos(arg) : wi=sin(arg) 

for i#=m# to n££t% step istep# 

j%=i%+nl% 

tr=wr*f r ( j%)-wi*f i ( 3%) : ti=wr*f i( ) +wi*f r ( j % ) 

fr(jS&) = fr(i*)-tr * fi(jsO-fi( i*)-ti 

fr ( i#)=fr ( £%)+ tr : f i ( ±%) =£ i ( ±%) + ti 

next 1% 

next 

nl%=istep% : goto 30320 
if itr%=1 then return 
for i#=1 to n££t% 

fr(i$g)*fr(ij6)/nfft$g : f i ( 1%) =f i ( i%) /nf f t% 
next i.% 
re turn 
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Signal or Sequence 


Fourier Representation 



(i ) 


(ii) 


(iii ) 


iX„ (jll) 


A 


- fl 0 o 

2 2 


ft 


Z(e J “) 


-101 




Mnterpreted as T (sec) 
(a) 


x a ( t ) 




°l- to —I 
Sampled 



i x(n) 



< 

► 
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Fig. 1 Fourier representations of (a) continuous- time signal and 
sequence derived from sampling, and of (b) finite-duration 
continuous-time signal, sequence derived from sampling and 
the sequence with implied periodicity. 
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Fig. 2 (a) Linear convolution of continuous- time signals and (b) 

circular convolution of sequences derived from sampling the 
signals . 
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Fig. 3 Schematic of experimental components with specimen showing the "stress 

wave factor" test configuration. (All tests in this study were conducted 
with transducers coupled face to face without a test specimen.) 
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Fig. 4 Schematic of digital signal processing procedures applied to experimental 
input and output sequences to evaluate unit-sample response and frequency 
response of experimental system. 
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Fig. 5 Time history, magnitude and phase of (a) input signal and 
(b) output signal. 
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Fig. 6 Experimentally determined (a) impulse response and (b) 
frequency response of test system. 
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Measured and predicted system responses: (a) 1.5 MHz center 

frequency tone burst input and corresponding (b) experimental 
output and (c) predicted output. 
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Fig. 7 
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Fig. 8 Measured and predicted system responses: (a) 2.5 MHz center 

frequency tone burst input and corresponding (b) experimental 
output and (c) predicted output. 
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Fig. 10 Measured and predicted system responses: (a) One 

cycle of 1.75 MHz frequency input at 2.3 Volts 
peak-to-peak and corresponding (b) experimental output 
and (c) predicted output. 
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